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Abstract 

Using both a recent calculation by Bruch of the damping of the motion of 
a monolayer nitrogen film oscillating harmonically on a metallic surface due to 
Ohm's law heating and a Thomas-Fermi approximation treatment of the Ohm's 
law heating mechanism, which accounts for the nonzero thickness of the surface 
region of a metal, it is argued that this mechanism for friction is able to account 
for recent measurements of the drop in the friction for a nitrogen film sliding 
over a lead substrate as it goes below its superconducting transition temperature. 
Bruch's calculation is also made more transparent by re-doing the calculation for 
a film sliding at constant speed, instead of oscillating. Using this treatment, it is 
easily shown that Bruch's calculation is equivalent to integrating Boyer's solution 
of the problem of a charge sliding over a metallic surface over the charge density 
of the monolayer nitrogen film. 

1. Introduction 
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Motivated by attempts [1-3] to explain a recent quartz crystal microbalance 
experiment [4] which shows a rapid drop in the friction of a film of nitrogen 
molecules sliding on a lead substrate, on dropping below the superconducting 
transition temperature Tg of the substrate, Bruch has recently done a calcula- 
tion of the electronic contribution to the friction for a monolayer nitrogen film 
executing simple harmonic motion on a metallic substrate [5]. In contrast to cal- 
culations of a single molecule moving on the substrate [6, 8], which require that 
the molecule possess a larger dipole moment or charge than is generally accepted 
for adsorbed molecules, in order to explain the experimental results of Ref. 4 
[1-3] , Bruch's results suggest that the field due to the quadrupole moment of the 
nitrogen molecule can explain the microbalance experiments [5] if the molecules 
form a monolayer film. The reason for this is that for a monolayer film the field 
inside the metallic substrate on which the film is moving falls exponentially to 
zero over a distance of the order of a lattice constant of the film below the surface 
of the metal. Since the distance over which the field is nonzero is much smaller 
than a mean free path, Bruch pointed out that the problem must be treated in 
the anomalous skin effect regime [7], in which only those electrons with velocities 
nearly parallel to the film remain in this region for a sufficient length of time to 
be significantly accelerated by the electric field. Consequently, only a fraction 
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{G£)~^ of the electrons (where G is the magnitude of a typical reciprocal lattice 
vector and £ is the mean free path) can be accelerated by the field. This results in 
an enhancement of the effective resistivity of the skin depth region, which leads 
to an enhancement of the rate of dissipation and the contribution to the kinetic 
friction due to Ohm's law heating. 

In addition to the Ohm's law heating mechanism for electronic friction put 
forward in Refs. 1 and 2 as a possible way of explaining the experimental results 
reported in Ref. 4, there is another mechanism for electronic friction (actually, the 
most commonly suggested mechanism for this phenomenon) which is due to the 
creation of electron- hole pairs of nonzero energy [9- 12]. The physical difference 
between these two mechanisms can be understood as follows: In the electron- 
hole pair mechanism for electronic friction, the energy loss due to friction is 
ascribed to the energy needed to create the electron-hole pairs. In contrast, in 
the Ohm's law heating mechanism considered in references 1, 2 and 5 and in 
this article, the energy loss need not be due to the energy cost necessary to 
create electron hole pairs. In this mechanism, if we consider the case in which 
the film is slid along at constant speed (e.g., by an applied force), the electric 
field resulting from the sliding film results in a screening charge near the surface 

of the metal, which is dragged along with the film. This results in an electric 
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current. Let us first consider only elastic scattering of the electrons (the dominant 
contribution to the resistivity well below the Debye temperature). During the 
sliding, electrons get scattered elastically by impurities and other defects. This 
would result in a reduction in the drift velocity, and hence the electric current, 
except that we force the current to remain constant by forcing the film to move 
at constant speed. In order to maintain the current, the electric field acting 
on the conduction electrons due to the film must accelerate them in order to 
maintain the drift velocity. The work done by this field is identified with the 
contribution to the dissipation produced by Ohm's law heating. In the Ohm's 
law heating mechanism considering only elastic scattering of the electrons, the 
electron-hole pairs resulting from the scattering of the electrons by defects in 
the substrate have zero energy because the scattering is elastic. In addition to 
the contribution due to elastic scattering, there is also inelastic scattering of the 
electrons by phonons, which results in an additional energy loss from energy 
transferred from the electrons to the phonons. Below Tg, the above scattering 
mechanisms do not occur for the superconducting electrons because of the gap in 
their excitation spectrum. The screening charge will be transported entirely by 
the superconducting electrons because they can fiow without electrical resistance, 
and hence are able to short circuit the current due to the normal electrons. 
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Since the superconducting electrons are not scattered, the Ohm's law heating 
mechanism for dissipation (and hence kinetic friction) described above does not 
operate. 

Persson has argued that the electron-hole pair mechanism should dominate 
over the Ohm's law heating mechanism by three orders of magnitude for a charged 
ion moving above the surface of a metallic substrate[13]. Since this mechanism 
depends on the density of normal electrons, which does not drop rapidly on 
falling below the superconducting transition temperature, however, it cannot 
account for the experimental result reported in Ref. 4. The calculation of the 
Ohm's law heating contribution of the friction, including the anomalous skin 
effect conductivity, as suggested by Bruch[5], however, can be of the same order 
of magnitude as the electron-hole contribution to the friction or greater. Since 
the resistivity drops to zero over a relatively small temperature range on dropping 
below Tc, the latter mechanism for electronic friction drops rapidly in much the 
same way as in Ref. 4, and, as we shall see, it is large enough to explain the 
experimental observations. 

In Bruch's treatment, the simple harmonic motion of the film results in an 
electric field with many harmonics of the frequency of the film's oscillations, which 
is an unnecessary complication in his treatment. In contrast, if one considers the 
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film to be sfiding over the substrate at constant speed, the field will possess 
only a single harmonic of the "washboard frequency," which is the characteristic 
frequency in this formulation of the problem. In section II, Bruch's calculation 
is reformulated by assuming that the film slides at constant velocity, rather than 
executing simple harmonic motion. It is expected that this will give the same 
value for the friction, and it simplifies the calculation. Furthermore, doing the 
calculation in this way makes it possible to verify that the force of friction acting 
on the sliding film calculated from the Ohmic heating inside the metallic substrate 
is equal to the friction found by calculating the force on the film due to its image 
field acting back on it, as is required by energy conservation. (This is not possible 
in Bruch's treatment because when averaged over the oscillation period of the 
film, the force of friction averages to zero.) It is argued in section III, that this 
version of Bruch's treatment is identical to Boyer's treatment [6], which is identical 
to the treatment of the problem due to Tomassone and Widom[8], used in Refs. 
1 and 2. Since Ref. 13 stresses the necessity of taking into account the nonzero 
thickness of the surface region, and since, Bruch's work is a classical treatment 
of the surface of a metal (i.e., one which assumes the surface region has zero 
thickness), in section IV, a calculation is presented which takes into account the 
fact that the surface region, when treated quantum mechanically, has a nonzero 
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thickness, in contrast to the zero thickness that it has in treatments of this 
problem using classical electrodynamics [1-3,5-8]. 

2. Bruch's Idea Applied to a Uniformly Sliding Film 

Let us take the region with 2; < to be occupied by the metal, with the 
region with z > occupied by free space. Bruch [5] writes the z-component of 
the electric field in terms of its Fourier transform on the x and y coordinates, 
E{G,z), 

£;,(r)= J]£;(G,2)e^G-, (1) 
G 

where G denotes a reciprocal lattice vector of the nitrogen film. Using the require- 
ment that the field have zero divergence (in regions with zero charge density), 
the field components parallel to the surface are given by 

E||(r) = ^(iG/G2)a£;(G,2)/a2e'«-'-. (2) 

G 

In regions in which there is no net charge density the current density J(r) can 
similarly be expressed in terms of the Fourier transform of its z-component. In 
Bruch's work, the film is assumed to execute simple harmonic motion of frequency 
as a rigid unit, which means that the time dependent fields are obtained by 
replacing r by r — Acos{Qt), where the vector A has magnitude equal to the 
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amplitude and direction in the direction of motion of the film. When this sub- 
stitution is made in Eqs. (1) and (2), the time dependence of the fields is a 
sum of harmonics Kfl, where K is an integer, with Bessel function coefficients 
Jk{G-A). In Bruch's solution it is necessary to deal with all of the time Fourier 
components. In this section, I propose that the inverse slip-time, which Bruch 
obtains, can be obtained much more easily by considering the damping of a film 
moving at constant velocity v instead. For the present treatment, in which the 
film slides at constant velocity v, we obtain the time dependence simply by re- 
placing r by r-vt, which results in a field which contain only a single time Fourier 
component for each reciprocal lattice vector, with frequency equal to the "wash- 
board frequency," G • v, and there is no sum over harmonic with Bessel function 
coefficients. Because of the linear nature of Ohm's law, this method must give 
the same value for the slip-time as Bruch obtains. Bruch obtains a relationship 
between the the time and space Fourier transforms of the z-component of the cur- 
rent density inside the film, and the z-component of the time and space Fourier 
component of the field on the surface of the metal, just inside the metal [Eq. 
(2.10) in Bruch's paper] by solving the linearized Boltzmann equation simultane- 
ously with Faraday's and Ampere's laws. Because the frequencies involved in this 
problem are quite small, Bruch solves these equations in the zero frequency limit. 
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As a consequence, each time Fourier component of the field and the current den- 
sity satisfies the same equations (since in the zero frequency hmit, the coefficients 
multiplying the fields in the equations are obviously independent of frequency). 
Thus, Bruch's solution of Boltzmann's equation with Faraday's and Ampere's 
law can be equally well applied to the present case of a film sliding at constant 
velocity, for which there is only one time Fourier component for each reciprocal 
lattice vector G, e.g., E(G, z) for the z-component of the electric field. One ob- 
tains for the relationship between the Fourier transforms of the z-component of 
the current density and field just below the surface of the metal 



where ug = Sa{l —p)/{4G£), where £ is the mean free path, p is the fraction 
of the conduction electrons which are specularly reflected at the surface of the 
metal at z=0 and a is the Ohm's law electrical conductivity. This is Bruch's Eq. 
(2.10). Substituting Eq, (3) in the standard boundary condition [14] 



J{G,z = 0-) = aGE{{G,z = 0-), 



(3) 



J,(r,2; = 0-) = (47r) 



-1 



d_ 
di 



[E,{r,z = 0-)-E,{r,z = 0+)], 



(4) 



we obtain 



Y,<^GEiG,z = 0-)e 



iG-r 



G 
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iA7r)-^d[E,{z = 0-) - E,iz = 0+)]/dt = 

J2{G ■ v)[E{G, z^O-)- E{G, z = 0+)]e^«-, (5) 
G 

where we have used the fact that E and J have the time dependence e~^^'^^ for 
a film sliding at constant speed, for the reasons given above. Eq. (5) can be 
written as 

(1 + iXG)E{G, 0") = E{G, 0+) = Bgs^' + Ei{G, 0+), (6) 

where -Ej(G, z) and E{G, z) are the Fourier transform of the z component of the 
contribution to the field from the film in the absence of the substrate and inside 
the metal, respectively, and B^e'^'^ is the contribution of the Fourier transform 
of the z component of the field outside of the substrate due to the charge density 
induced by the film. Bq is a constant to be determined by the boundary con- 
ditions. The parameter Xq = A7raQ/{G ■ v), where (G • v) is the "washboard" 
frequency of the film. In addition to Eq. (6), we must require continuity of the 
component of the field parallel to the surface, which is given by 

dE{G,z)/dz\,=o- = [d{BGe^')/dz + dEi{G,z)/dz]\,=o+, (7) 

which using the result 

E{G,z) ^ E{G,z^O-)e^' 
11 



given in Eq. (A13) of the appendix of Ref. 5 (It follows from the solution of 
Ampere's and Faraday's law in the zero frequency limit.) gives 

-E{G,0+)=Bg-E,{G,0-). (8) 

The solution to Eqs. (6) and (8) is 

E{G, 0+) = [2/(2 + iXG)]Ei{G, 0"), (9a) 

and 

Bg = [«Ag/(2 + tXG)]Ei{G, 0") ^ (1 + 2t/\G + ...)Ei{G, 0"), (96) 
where Ei{G, Z) can be found using 

Ei{r) = J d'r'p{r)j^^, (10) 

where p(r) is the film's charge density. It can be crudely modeled by three charges 
along the axis of a molecule ip, in the small £p limit, two charges of charge +q 
at the outer edges of the molecule and a charge -2q at its center. The value of q 
is chosen such that the resulting quadrupole moment has the experimental value 
9. Then, in the unit cell near the origin 

p(r) = q5{z + h)J2 [^^'^ (r - P/3 - ^p) + S^^^ {r - + ip) 

(3=1,2 
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-25(2)(r-p^)], (11) 

where 5^^) denotes a two dimensional delta function and pp denotes the position 
of a molecule in the unit cell. Writing 



r — r 



r — r 



and substituting in Eq. (10), we obtain 

Ei{G,0-) = -(7r^Me)e-^'^5^e-^G-^'^(G-i;3)2, (12) 

where pp is the location of the molecule in the nitrogen film unit cell, £j3 is 
the symmetry axis of molecule, and 9 is the quadrupole moment of a single 
molecule. The inverse slip-time found by calculating the force exerted on the film 
by the image field outside the metal, 

J2iiG/G^)d[BGe^'e'^-'']/dz, 
G 

is given by (NMv)~^ times the force or 

(NMvy^ I d\p{r)Y^{iG/G'^)d[Be^'e''^-'']/dz. (13) 

G 

From Eqs. (9b) and (13), we obtain for the inverse slip-time 



= -T-. )iy e-2^'^(4G£)(G-v)(G/G') 

G 
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\Y,{G-ipfe'''-P^\\ (14) 

which gives 

r-i = I ) y (4G'£)(G'2/G')e-2G/^ 

^3(1 -p)aAcM^ ^' ^ 

G 

I J](G-4)V«-''^|2, (15) 

where we have taken v to be along the x-axis. By energy conservation, this 
expression for r^^ must be equal to the value calculated by setting the power 
loss due to Ohm's law heating equal to NMv^/t^ the rate at which the viscous 
electronic contribution to the friction force, NMv/r, does work on the film. 
Using the method used in this paper of assuming that the film is sliding at a 
speed V, rather than oscillating (as assumed in Bruch's paper), the Ohm's law 
heating power is given by 

P = 2NAcJ2G~'\Ei{G, 0-)|2(G ■ v)/ctg (16) 

G 

The factor NAc appears because when P is calculated by integrating J(r) ■ E(r) 
over the volume of the metal, the volume integral involves evaluating the integral 

G G'. We obtain with this proceedure[15], substituting Ei 
from Eq. (12) in Eq. (16), the value for obtained in Eq. (15). 
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3. Treatment of Electronic Friction using Boyer's Solution 

Boyer[6] solves the problem of a charge or electric dipole moving above (i.e., 
outside the metal) and parallel to the surface of the metal. He solves the electro- 
dynamics problem subject to the same boundary conditions as Bruch uses [14], 
[Eq. (4) above] , which for Boyer's problem is written as 

J^ = aE= {ATT)-^dE^/dt = {Att)-\ ■ VE^. (17) 

When terms of up to first order in v are kept, we obtain Boyer's result for the 
electric field. The force exerted by this field on the moving charge gives the force 
of friction due to Ohm's law heating in the metal. The force of friction found 
by the formalism due to Tomassone and Widom[8] gives the same friction and 
hence is believed to be equivalent to Boyer's calculation. 

Boyer[6] finds, in addition to the electrostatic field, a contribution to the 
electric field above the substrate linear in the velocity v of a point charge q 
sliding above the substrate, given by 

E(r) = -(W27ra)^[^f^], (18) 

assuming that the sliding velocity is in the x-direction, where r' = vtx + hz is 
the location of the moving charge (where x and z are unit vectors in the x and 
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z directions, respectively). In order to apply this result to a monolayer film of 
charge density p(r), let us multiply Eq. (18) by p(r) and integrate over volume 
to obtain E(r), 



where we have written the field in Eq. (19) in terms of its Fourier transform. For 
a periodic monolayer film a height h above the surface of the substrate, p(r) has 
the form 



where G denote the reciprocal lattice vectors of the film and pc = -^c^ d'^fP2{r)e 
where the u on the integral sign signifies an integral over a unit cell of the film, Ac 
is the unit cell area and p2(i') is the charge per unit area of the film. For simplic- 
ity, we are modeling the film by a collection of point charges. If the nonzero size 
of the charges in the film were taken into account, there would be form factors 
introduced in the summations over G, which would fall off rapidly with increas- 
ing magnitude of G. This can be approximately accounted for by including only 





(19) 




(20) 



G 
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the first one or two terms in the sums. Substituting Eq. (20) in Eq. (19) gives 
E = 2v{27:a)-^Y.PG J dk,{j^^)e^-^'n-vtx-{z+h)z]^ ^21) 

G ^ 

where ry is the projection of r in the plane of the film, and where the x and y 
components of k are equal to the x and y components of film reciprocal lattice 
vectors. For components of E parallel to the substrate we obtain on performing 
the integral over 

E(r) = {v/a)J2PG{GGa:/G)e-^^'+''^e'^^-''+^^''^\ (22) 
G 

To find we model the charge density of of each molecule by three charges 
as was done in the last section. Then, substituting p(r) from Eq. (11) in the 
integral for pc under Eq. (20), we obtain 

PG = -4gA-i Yl ^^^'(G ■ £^/2)e^G-P^. (23) 



Therefore, 

E(r) = -v4qA-^a-^ ^^^(G • 1/2) 

AG 

X {GGa,/G)e-^^'+^^e'^^<P^+''^-^-''*^ , (24) 

which reduces to an expression resembling Bruch's for the field in the small £p 
limit with the quadrupole moment of the molecule 6 equal to g£^, if we pretend 
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that the system is not in the anomalous skin effect regime, and replace uq by cr 
for the purpose of making a comparison. The force of friction acting on the film 
is given by 

F = (l/2)i?e J dVp*(r)E(r) = -{NO^/a) 
x^c"' E I ■ ^>)2e^«-^^|2(GG,/G), (25) 

G 

where N is the number of molecules in the film, and thus the inverse slip-time 
is given by 

= (F/NMv) = 

A-\e^/Ma) J2 I • 4)'e^«•^''|2(G'2/G')e-2«^ (26) 

G 

where M is the mass of an adsorbed molecule and £^ is a unit vector along £f^. If as 
discussed above, we replace ctq of section II by u, Eq. (15) becomes identical to 
Eq. (26). This demonstrates the equivalence of Bruch's and Boyer's treatments. 

4. A Treatment of the Problem which Includes the Nonzero Width of 
the Surface Region 

Both Bruch's treatment of the problem and Boyer's are based on the classical 
model for a metallic surface, in which the electronic charge density drops to 
zero immediately on leaving the metal. In a quantum mechanical treatment, in 
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contrast, the electronic charge density drops to zero over a distance of the order of 
a couple of Angstroms. The film almost certainly resides in a region just outside 
the bulk of the metal, in which the electronic charge density is nonzero (although 
it is decaying exponentially here). Ying, et. al., studied electron screening in 
the surface region using the Thomas-Fermi approximation [16]. They find that a 
charge placed at the surface will be screened, with a screening length not much 
longer than that in the bulk metal. The bulk metal Thomas Fermi expression for 
the screening charge density of a point charge can be used to model the charge 
density in the surface region analytically if one fits the screening length to that 
found in Ref. 16. The bulk Thomas expression for the screening charge density 
for a point charge is 



where q is the point charge whose screening is being considered and ks is the 
inverse Thomas Fermi screening length, which we will take here to be a parameter 
to be fit to Ying, et. al.'s calculations [16]. In Ref. 16, the quantity 



p,(r) = -(47r)-ig/c2e-^^7r 




(27) 




(28) 
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is plotted as a function of z. The model of Eq. (27) gives for this quantity 

-(l/2)g/c,e-^^l^l, (29) 

which strongly resembles the quantity plotted in Fig. 1 of Ref. 16. Then to fit the 
present approximate model to the results of Ref. 16, we can simply choose a value 
of ks for which Eq. (26) reproduces each of the plots in Fig. 1 of that reference. 
This method is rigorous when one can use the quasi-classical approximation, 
which is accurate if typical values of the screening length are much smaller than 
the thickness of the surface region. (See the appendix.) When this is not an 
appropriate limit, for example if the film is further out from the bulk of the 
metal, there is asymmetry in the screening charge density [17], but it is not such 
an extreme asymmetry as to invalidate using a spherically symmetric screening 
charge density, as is done here as a first approximation. 

On the basis of this model, we conclude that when a film with charge density 
p(r) moves along the surface of the metal with a velocity v, the screening charge 
density given by 

p,(r) = -{2n)-'kl J dV J d'k^^^pir') (30) 

moves with the same velocity resulting in a current density Jg = vps{r). Substi- 
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tuting for the p(r') in terms of its Fourier transform, we obtain 

G * 



(l/2)fc,vA, (G2 + ;,2)l/2 (31) 

The resulting Ohmic heating contribution to the force of friction can be found 
from 

-1 / IT |2j3 



FfricV = a-' J iJsl'd-'r. (32) 

Substituting Eq. (31) into Eq. (32) while replacing a by an effective conduc- 
tivity (Jg ~ cr/G£ to account for the anomalous skin effect [7] (like the effective 
conductivity used in section II) and placing it within the summation over G, we 
obtain for the force of friction 



Ffric = AkjvJ^^G'iG^ + k!r^^^\PG\^- (33) 

G 

Then, we have for the inverse slip-time 

—1 -Ffric 
T = — 



NMv 



24 M^"G^^G2 + /.2)-3/2|^^G.i,)VG-^^P, (34) 
G f3 

where we have substituted for pQ using Eq. (23), where A is the area of the film 
Ac is the area of a unit cell and N is the number of unit cells in the film. For 
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kg ~ 10^cm~^ and 9 ^ 10~^^esu, we obtain 10^^ s~^, which is much larger 

than the electron-hole creation mechanism discussed earlier [13]. The quantity 
T^^ found from Eq. (34) is about four orders of magnitude larger than r^^ found 
from Eq. (15) because Eq. (34) does not contain the factor e~'^^^, which appears 
in Eq. (22) (where G = 2a~^ and h = 2A). If we use a value for ks a factor of 5 
smaller (which is not unreasonable considering that we are in the surface region 
where kg could be noticeably smaller than its bulk value because the conduction 
electron density is smaller here), r"^ becomes of order 10^s~^, comparable to the 
experimental value[4]. 

Bruch's use of a classical (i.e., zero thickness) surface would be valid if the film 
resided well above the surface where the (classical) method of electrical images 
should be valid[17]. The treatment in this section assumes that the film resides in 
a part of the surface region in which the electron density is closer to its bulk value. 
The two treatments bracket the true situation, in which the film lies between 
these two extremes. Since both treatments give a large enough magnitude for 

to account for the experimental results, one can say with confidence that the 
Ohm's law heating contribution to the friction is of sufficiently large magnitude 
to account for the experimental results of Ref. 4, as put forward in Refs. 1 and 
2. 
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5. Conclusions 

Bruch's calculation of the electronic friction acting on a film of nitrogen 
molecules harmonically oscillating on a metallic substrate, is re-done for the 
simpler case in which the film is sliding at constant speed. This treatment, 
which should give the same value for the slip-time as Bruch's treatment, clarifies 
Bruch's treatment and allows one to easily demonstrate that the force of friction 
found by calculating the force of the image charge acting back on the film and by 
calculating the Ohm's law heating inside the metallic substrate are equal. It also 
allows one to demonstrate the equivalence of Bruch's treatment with that due 
to Boyer of a charge sliding over a metallic surface. Since Boyer's treatment is 
equivalent to that used in Refs. 1 and 2, one is confident in saying that Bruch's 
treatment of the problem is an extension of these methods, which allows one to 
include the anomalous skin effect. 

Persson has argued that the contribution to the friction due to the creation 
of electron-hole pairs should dominate over the Ohm's law heating mechanism by 
three orders of magnitude for a charged ion moving above the surface of a metallic 
substrate[13]. Since this mechanism depends on the density of normal electrons, 
which does not drop rapidly on falling below the superconducting transition 
temperature, it cannot account for the experimental result reported in Ref. 4. 
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The calculation of the Ohm's law heating contribution of the friction presented 
in the previous sections, including the anomalous skin effect conductivity, as 
suggested by Bruch[5], however, can be of the same order of magnitude as the 
electron-hole contribution to the friction. Since the resistivity drops to zero over 
a relatively small temperature range on dropping below Tg, the latter mechanism 
for electronic friction drops rapidly in much the same way as in Ref. 4, and, as 
we have seen, it is large enough to explain the experimental observations. 
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Appendix: Quasiclassical Treatment of the Screening of a Charge in 
the Surface Region 

In the surface region, where the electron charge density is decreasing from its 
bulk value down to zero, the wave functions in the jelleum model must take the 
form 

V'k(r)=e*-7k(^), (^1) 
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where /k(-z) drops from the value it has for z in the bulk region to zero for z well 
above the surface. Then the electron charge density is equal to 

E^k|V'k(r)|2 = 5^nk|/k(;^)|2, {A2) 

k k 

where is the Fermi function, {e^^^^)~^^)I^BT _|_ where e(k) is the electron 
energy of wavevector k and is the chemical potential. The potential 0(r) 
satisfies Poisson's equation 

VV(r) = -47r5^nk|/k(^)p. (^3) 
k 

Let us now consider the screening of a point charge q placed at a point r = zqz 
inside the surface region, where 5 is a unit vector in the z-direction. In the 
linearized Thomas-Fermi treatment of screening, we assume that (j){r) is changed 
by a small amount 5(j){v) because of the point charge, add this change in the 
potential to jj. and linearize in 5(j). Carrying this out we find that 5(j){v) satisfies 

V25(/)(r) = -kl{z)5(t){Y) + 47r5(r), (A4) 

where 

/c2(^)«47rj]5(e(k)-/.)|/k(z)|2 (^5) 

k 

in the low temperature limit. Taking the Fourier transform of Eq. (A4) with 
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respect to the components of r parallel to the surface it becomes 

(f5(f){^,z)/dz^ = -[kl[z) - /c^)5(/)(k||, ^) + Att5{z - zq), (A6) 

where (^0(k||, z) is the Fourier transform on 5(j){r) over the coordinates parallel to 
the surface. Let us attempt to find a solution to Eq. (A6) of the form 5(j) = e~^^^\ 
Substituting in Eq. (A6) we obtain 

[{dS/dzf - d^Sjdz^ - {kl{^z) - /cpe-^ = 47r5(^ - z^). (yl7) 

If typical values of the screening length are small compared to the width of the 
surface region, we can for most z neglect d'^S/dz^ compared to (dS/dz)^. With 
this approximation, we obtain a solution 

z 

S(l) = {4:nq/2k',{zo)exp[-\ J k',{k')dz'\], {AS) 

where (kg)"^ = kg—k'^^. In the extreme limit in which typical k'g~^ values are much 
smaller than the width of the surface region, we can replace the integral in the 
exponent by kg{zQ){z — zq) to a good approximation. The resulting form for S(j) 
gives an inverse Fourier transform on k^^ proportional to 

-ks{zo)\r-zoz\ 

iA9) 



\r - zoz\ 

the form of the result for Thomas-Fermi of a point charge in the bulk of the 
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metal. While the extreme limit of screening length much smaller than the width 
of the surface region is not likely to occur, we do not expect the screening in real 
situations to be so qualitatively different than the bulk Thomas Fermi screening. 
For example, one would not expect the screening that occurs parallel to the 
surface to be qualitatively different than that normal to the surface, although 
there is certainly likely to be some anisotropy[17]. 
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